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Abstract 



We present a comprehensive analysis of a linear growth model, which com- 
bines the characteristic features of the Edwards-Wilkinson and noisy Mullins 
equations. This model can be derived from microscopies and it describes the 
relaxation and growth of surfaces under conditions where the nonlinearities 
can be neglected. We calculate in detail the surface width and various cor- 
relation functions characterizing the model. In particular, we study the 
crossover scaling of these functions between the two limits described by the 
combined equation. Also, we study the effect of colored and conserved noise 
on the growth exponents, and the effect of different initial conditions. The 
contribution of a rough substrate to the surface width is shown to decay 
universally as Wi 

(0)(6/e(0)'^/', where ^ ~ t^/' is the time-dependent 
correlation length associated with the growth process, Wi{0) is the initial 
roughness and the correlation length of the substrate roughness, and d 
is the surface dimensionality. As a second application, we compute the 
large distance asymptotics of the height correlation function and show that 
it differs qualitatively from the functional forms commonly used in the in- 
tepretation of scattering experiments. 
PACS numbers: 68.55.-a, 68.10.Jy, 05.40.+j. 



1 Introduction 



The dynamics of interfaces ranging from dendritic growth to flame front 
propagation can often be described by relatively simple evolution equations 
jT] . Typically such evolution equations are given in terms of partial differen- 
tial equations with a stochastic noise component. Perhaps the best known 
example is the nonlinear Kardar-Parisi-Zhang (KPZ) equation which 
describes kinetic roughening of randomly driven interfaces such as growing 
surfaces or flame fronts in forest fires Q . Complete understanding of these 
nonlinear equations is still mostly lacking. 

Under certain circumstances discussed below the relevant nonlinearities may 
be so weak that a fully linear model can provide an adequate description. 
One example is the growth model of Edwards and Wilkinson (EW) Q, 
which describes the sedimentation of granular particles under gravitation. 
Another important model is the noisy Mullins equation discussed by Wolf 
and Villain and others [|, ||, |7) (the Mullins-Wolf- Villain (MWV) equation) 
in the context of molecular beam epitaxy (MBE). Being linear, both of 
these equations have been analyzed in some detail. However, recently it 
has been shown both from macroscopic arguments and through more 
microscopic derivations ^ that for some cases involving surface diffusion 
and desorption, a more general linear equation of the form 

dth = uiV^h + U2V^h + r/ (1) 

emerges, where h = h{x, t) is the surface height above a d-dimensional 
substrate, rj{x, t) is a noise term, and vi and 1^2 are parameters. Since the 
gradient terms of Eq. ([l|) are simply a combination of the EW and MWV 
equations, we call it the Combined Linear Growth (CLG) equation. Stability 
requires that vi > and 1^2 < 0. While it is physically possible that z/i < 
(see Sec. 2.2), the treatment of this case requires the inclusion of additional 
nonlinear terms in ([l|), and will not be addressed here. 
The purpose of the current work is to present a detailed analysis of Eq. 
(1), which is missing so far. This is useful for two main reasons. First, the 
calculations in this work generalize the previous results obtained for the EW 
and MWV equations which are somewhat incomplete and scattered in the 
hterature g |, |, |, |ll|, ||, |l|, Second, Eq. ^ is the simplest 



example of a growth equation with an intrinsic length scale. Balancing the 
two gradient terms in (||) one finds that they become comparable at the 
scale 

r = v/|^2lM. (2) 
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The kinetic roughening process is governed by the fourth derivative term on 
scales smaller than i* but the second order term dominates on scales larger 
than i* . Physically, the two terms represent different relaxation mechanisms, 
through surface diffusion (fourth derivative), evaporation-condensation or 



step edge barriers (second derivative) 17|. Of course, writing a continuum 
equation with an intrinsic scale is meaningful only if this scale much exceeds 
the microscopic cutoff, given by the lattice spacing a; the detailed estimated 
of ui and V2 derived in Sec. |2| show that i* is indeed mesoscopic, in the 
sense that i* ^ a, under typical conditions. Mesoscopic intrinsic length 
scales such as domain or terrace [19| sizes play an important role in the 
kinetic roughening of real surfaces. Moreover, the competition and interplay 
between different relaxation and roughening mechanisms is probably typical 
in many experimental situations; for this reason combined linear equations 
like (|^) have already been used extensively in the analysis of experimental 
data @. 

A major advantage of working with linear growth equations is that they 
allow us to explicitly compute any statistical quantity of interest, rather 
than just extracting the values of scaling exponents, which have been the 
focus of most previous studies of kinetic roughening |21|. We will exploit 
this fact to address two questions of direct experimental relevance: The 
evolution of the initial substrate roughness during growth, and the shape 
of the height-height correlation function. In both cases we find that the 
heuristic expressions commonly employed in the experimental literature are 
incompatible with the explicit calculations. 

The organization of the paper is as follows. We shall first briefly review 
the physical background and various arguments leading to Eq. ([l|), and in 
particular various interpretations of the coefficients vi and 1/2- Following 
this. Sec. ^ contains a full solution of the CLG model in terms of the 
relevant measures of the surface roughness. We concentrate on the scaling 
behavior and finite size dependence of the surface width w{L,t), the equal 
time height-height correlation function G{x,t), and the saturated height- 
height correlation function Cs{x,t). In particular, we study in detail the 
crossover behavior of these quantities to the well known limits given by the 
EW and MWV equations. In Sec. 4 we discuss the experimentally important 
effects of substrate roughness, the shape of the height correlation function 
and the influence of long ranged noise correlations. Finally, summary and 
conclusions are given in Sec. 5. 



2 Derivation of the Linear Growth Model 

In this section we provide some microscopic justification for the combined 
linear growth equation ([l|). We consider two different physical situations, 
corresponding to a surface in thermal equilibrium (Sec. |2.1[ ) and a vicinal 
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surface growing in the step-flow mode (Sec. 2.2), respectively. In both cases 
the fourth order derivative term in reflects capillarity-driven surface 
diff'usion [jl^ , while the second order term will be seen to arise from distinct 
mechanisms. 



2.1 Equilibrium Dynamics of a Solid— on— Solid Model 

We consider a one-dimensional solid-on-solid surface described by a set of 
integer height variables hi defined on a lattice. The energy of the surface is 
given by the Hamiltonian ||2^ 

n = {K/2)J2\h^-h^+l\. (3) 

i 

The surface evolves according to the following dynamic processes Q: Par- 
ticles are deposited {hi — > /ij + 1) at a constant rate F; they evaporate 
{hi hi — 1) at rate 

Wr = ko exp[- {E's + niE'j^)/kBT] (4) 

and jump from site i to site j = i ±1 at rate 

= ko expi-iEs + n^EN)/kBT]. (5) 

Here /cq defines some microscopic hopping rate of the order of a typical 
phonon frequency. We will generally measure time in units of k^^ , so that 
effectively k^ = 1; likewise the basic length unit will be provided by the 
lattice constant. The energy barriers in (^) and each have a substrate 
contribution {E'g and Es) and a bonding contribution proportional to the 
lateral coordination number 

Ui = 9{hi+i - hi) + 9{hi^i - hi) (6) 

which takes the values rii = 0, 1, 2 in one dimension. Detailed balance rela- 
tive to (^) holds if 

En = E'n = K and F = ko exp[-{E's + K)/kBT]. (7) 

The distinguishing feature of the diffusion rates (|5|) is that they depend only 
on the environment at the initial site i. It was shown elsewhere |2^ how 
this fact can be used to exactly derive the continuum equation of motion 
for the surface, in the case where only surface diffusion is allowed. Here we 
generalize the approach of ||2^ to include desorption and deposition. Note, 
however, that this derivation is only valid when the surface is in equilibrium 
with the vapor, as expressed by the second condition in (0). 
From the master equation of the process one easily derives the following 
equation of motion for the ensemble averag ed height g, |2^, 

±{h,) = ie-^«/^-^(v2A). - e-^'s/'-^X. + F (8) 
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where (V^A)i = Aj+i + Aj_i — 2Xi denotes the lattice Laplacian, and Aj = 
(exp[— (iir/fcBT)nj]). This quantity is related to the local chemical potential 

^ii 11, 

\i=eM-{K-^ii)/kBT]. (9) 

We now pass to the continuum limit (hi), fii h{x, t), fj,{x, t), where h{x, t) 
and fJ-{x,t) are averages taken over some large region (still small on the 
macroscopic scale) centered around i = x. The local chemical potential 
^{x, t) is then determined by the local surface curvature via a Gibbs-Thomson 
relation 

^ = -7(V/i)V2/i (10) 

where the stiffness 7 is a nonlinear function of the local surface slope that can 
be directly computed from the Hamiltonian (0) p3|; at zero tilt (V/i = 0) 



7(0) = kBT[cosh{K/2kBT) - 1] . (11) 

For slowly varying, macroscopic profiles the typical curvatures are small, so 
that d) can be expanded in /i. This results in the macroscopic equation 

dth = -^{kBT)'^e-^^S+K)/kBTy2^^(yf^^y2^^^ (^2) 

+(A:Br)-ie-(^s+^)AsT^(V/j)v2/i 
which can be used e.g. to predict the decay of periodic surface modulations 



1 23]. Note that the deposition term has disappeared due to the second of 
the detailed balance conditions (0). Because of the orientation dependence 
of the stiffness, ( |l2|) is highly nonlinear. 

Here we are primarily interested in the mesoscopic fluctuations around a sur- 
face that is on average flat. Eq. (|12|) can then be linearized by expanding 7 
around the average orientation u = ( V/i) . It is not strictly speaking consis- 
tent to keep only the terms linear in h, since the nonlinearities arising from 
the expansion of the stiffness in the second term of (|T^ ) (such as (V/i)^V^/i) 
are more relevant in the renormalization group sense than the linear fourth 
order term; however we will ignore this difficulty in the interest of obtaining 
an analytically tractable model. Thus we arrive at the two systematic terms 
in dl]), and identify the coefficients as 

= {kBT)-mu)e-(^'s-^''y'-^, U2 = -l(A:Br)-i7(^)e-(^^+^)/'=-^. 

(13) 

To complete the derivation of (|l]), the statistics of the noise term has to be 
specified. This requires no further information, since detailed balance forces 
the stationary distribution of the continuum field h{x, t) to be governed by 
the Hamiltonian 

Tic = {l{u)/2) I dx {Vhf, (14) 
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as can be seen from a central limit argument applied to the sum of in- 
dependent local slope variables (^). A straightforward way to ensure the 
stationarity of exp[— "Hc/^s^] is to modify the distribution functional of the 
white noise into the following form: 

where 5 = -{^{u)/2kBT) f dk J dt{i^i+i^2k^)-^\fi{k,t)\^, and Z = J Vr] e^^'^\ 
where fi{k,t) is the Fourier transform of r]{x,t). This leads to the noise co- 
variance 

{7]{x, t)rj{x', t')) = 2{kBT/^{u)){ui + U2V^)5{x - x')5{t - t'). (16) 
Our derivation provides a microscopic basis for the classical theory of Mullins 



1 17], who showed that second and fourth order derivatives of h arise from 
evaporation-condensation dynamics and surface diffusion dynamics, respec- 
tively. We may further conclude that the length scale beyond which the 
lower order derivative associated with evaporation-condensation dynamics 
dominates is given by 

t ~ ^-V2/vi ~ exp[(l/2)(S;, - Es)/kBT] (17) 

in units of the lattice constant. Since typically the activation energies for 
evaporation much exceed those for surface diffusion, i* can be quite large 
at moderate temperatures. It should be noted, however, that due to de- 
tailed balance the length scale (|T7|) cannot appear in any stationary, equal 
time correlation functions, since these only depend on the Hamiltonian (|l^). 
Nevertheless, £* and a related time scale will appear in the time-dependent 
quantities, to be discussed in subsequent sections of this paper. 



2.2 Step-Flow Growth on Vicinal Surfaces 

Technologically important deposition techniques such as molecular beam 
epitaxy (MBE) are typically carried out at temperatures where desorption 
is negligible, so that effectively i^i = in (p^). However, as was first pointed 
out by Villain [^, under growth conditions other mechanisms related to 
growth-induced surface currents exist which generically give rise to a 
second order derivative in the continuum equation. A remarkable feature 
of such currents is that they can be destabilizing, leading to z^i < in (||). 
In the present work we focus on the kinetic roughening of a stable surface 
with vi > 0, and therefore we describe here only the simplest microscopic 
mechanism for the generation of a positive vi term through the "Schwoebel" 
effect involving step edge barriers [^] (for some other mechanisms see [^] ) . 
As in the preceding section we restrict ourselves to a one dimensional surface. 
In two dimensions the mechanism described here gives rise to an anisotropic 
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Laplacian z^||9y + i^±d'j_ in (||), with different coefficients z^y and parallel 
and perpendicular to the surface steps, at least one of which is negative p6| . 
We consider a vicinal surface with uniform step spacing i, which is assumed 
to be much smaller than the diffusion length i o governing the island spacing 
on a singular surface [^] ; this ensures that island nucleation on the terraces 
can be neglected, and the surface grows in the step flow mode ||2^ . Moreover 
we assume strong step edge barriers which effectively suppress any interlayer 
transport. Under such conditions every atom that is deposited on a terrace 
attaches to the ascending step edge, and the surface current is simply J = 



F(./2 (as before, F denotes the deposition rate) 24, 26]. The coefficient 



vi is then given by the negative derivative of J with respect to the surface 



inclination l/£ [24]. This yields 

vi = (18) 

Provided the capillarity-driven surface diffusion is not too strongly affected 
by the deposition, the expression ( p^ ) for 1^2 is still expected to be valid. 
Thus the crossover length scale (0) can be estimated as 



t^l-^^-v2/F^^U^ (19) 

where ^cap ~ {—V2/FY/^ is a length scale gauging the relative importance 
of capillarity and deposition ]p6[] . For an order of magnitude estimate, we 
note that 1^2 can be directly measured from the decay time of periodically 
modulated surface profiles. For semiconductor surfaces a typical value is 
|27] —U2 ~ l(/im)'* per hour, implying that ^cap ~ 1000 Aat a deposition 



rate of 1 s~^. This much exceeds the step spacing on typical vicinal surfaces, 
and thus, as in the case of equilibrium dynamics (Sec. |2.1| ), there are good 
reasons to expect i* to be large compared to the lattice constant. 
Of course, the most prominent effect of deposition is to provide an additional 
source of "shot noise" fluctuations. Assuming, again, that the (volume con- 
serving) fluctuations due to surface diffusion are not much changed by the 
deposition flux, we obtain the noise covariance 

(r/(x, t)ri{x',t')) = [F + 2{kBTv2h)y'^]Kx - x')5{t - t'). (20) 

A comparison between the strengths of the two components of the noise, i.e. 
F and kBTv2/l^ defines a further crossover length scale 

r* ~ {D,,xx/Ff'^ (21) 

where -DcoU = ^oexpj— (Sg + K)/kBT] defines, within the present model. 



the collective surface diffusion coefficient J23|. In contrast to the detailed 



balance situation of Sec. 2.1, here the two crossover lengths i* and £** need 



not be equal. Their ratio is of the order 

r/t* ~ {kBTyH-^j^/^ ~ exp{K/AkBT) (22) 
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where in the last step we have used the expression (^) for small T. Thus, at 
low temperatures I* ^ This provides some justification for neglecting 
the conserved noise component in (pO|), as will be done throughout Sec. |3[ 



3 Solutions for Various Physical Quantities 

In this section we summarize our results for the physically interesting mea- 
sures of the surface roughness of the CLG model, for arbitrary surface di- 
mensionalities d < z, where z is the dynamic exponent. The physically 
interesting quantities that we calculate are the surface width w{L, t) and 
two correlation functions G{x,t) and Cs{x,t). The surface width is the size 
of typical height fluctuations around the mean h = {{{h{x , t)) x) ri- 

w\L,t) = {{{h{x,t)-hf)x)r,, (23) 

where {■)x and (•)^ denote averaging over space and noise, respectively, and L 
is the lateral extent of the surface, assuming periodic boundary conditions. 
The two correlation functions can be derived from the general two point 
correlation function Cg{x,t,t'), which is defined as: 

Cg{x, t, t') = {{{h{x + x',t + t') - h{x', t')f)x')^ = Cg{x, t, t') . (24) 

All the correlation functions appearing in this section are thus dependent 
only on the magnitude x = \x\. The equal time correlation function G{x, t') = 
Cg{x, t = 0, t') and the saturated correlation function Cs{x, t) = Cg{x, t, t' S> 
ts), where tg is the saturation time to be defined later. We point out that 
the translational invariance of Eq. (||) makes the averaging over noise and 
space interchangeable when the noise correlations are also translationally 
invariant. In this section we assume Gaussian white noise of the form 

{r]{x,t)r]{x',t'))r^ = 2D5'^{x-x')5{t-t') (25) 
(r?(f,t)), = 0. (26) 

Later on in Sec. 4 we discuss the infiuence of the conserved noise compo- 
nent that appears in Eqs. (|lD and (|0|). Here we also assume a flat initial 



condition, h{x, t = 0) = 0. Growth on an initially rough surface is dealt 
with in Sec. |^. 

The results for the CLG equation trivially reduce to the limits of the EW 
or MWV equations when 1^2 = or vi = 0, respectively. Also, simple power 
counting shows that for an infinite system, the EW behavior will always 
eventually dominate. In this work, however, we are interested in investigat- 
ing the crossover time tc from MWV to EW growth. The crossover always 
occurs for an infinitely large system where the MWV behavior is dominant 
for early times (t <C tc)- Moreover, in a finite system for a suitable choice of 
the crossover length scale £*, the MWV growth can be made dominant for 



7 



all times. As was noted already, we focus on the stable case vi > 0, 1^2 < 0. 
If vi < 0, the early time growth can be of MWV type but the long time 
behavior of the surface is unstable. 

The main results for the CLG growth equation are given in the subsequent 



paragraphs. As usual |21], the asymptotics of the surface correlations involve 
the roughness exponent the dynamic exponent z and the exponent ratio 
(3 = xl z which describes how the surface width increases with time. The 
two limiting cases of ([l|) are characterized by the exponents 

zi=2, /3i = (2-d)/4, xi = (2-(i)/2 (EW) (27) 

as was first derived in |ll| , and 

^2 = 4, /32 = (4 - (i)/8, X2 = (4-d)/2 (MWV), (28) 

compare to |5|, ^, |l2| . It is understood that /? = x = implies logarithmic 
roughening. 

The crossover time scale is given by tc = \v2\/i'i and the saturation time 
(for finite system size L) by tg = L^'^ j {L^^v\ + |z^2|)- The crossover length 
scale t = \/\i'2\/i^i was defined in Sec. 1. In the following we use the 
dimensionless scaling variables: pi = vit/L"^, p2 = \i'2\t/L^, yi = vxtjx^ 
and 2/2 = 11^21^/2;^. The dynamic correlation lengths for the EW and MWV 
cases are defined as ^1 = (2z^it)^/^ and ^2 = (2|z^2|i)^''^- 

3.1 Surface Width 

The scaling function for the surface width is obtained from Eq. (p3|): 

on on 

w\L,t) = — l2xifJ,(pi,P2) = ^L'^'F^{p,,p2) , (29) 

where 

J'^m,P2) 2(2^)'^ A ^ ^ ^pik^+p2k^' ^ ' 

where i^d-i is the surface area of a d dimensional sphere, and j = 1, 2. We 
obtain the following power law behavior in the different time regimes: 

If pi » P2, then 

' L^xi , for f > ; 
t^/^i , for < t < . 



If pi <^ p2 then 



^2/^2 , for t < tc • (31) 




for t ^ ts ; 
for t<^ts. 



The main result here is that for a finite system the MWV behavior can 
dominate for all times, including the saturated regime. This is shown in 
Fig. 1, where it can be seen that even for pi = P2, the EW region can be 
made to vanish by choosing tg = tc- 
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3.2 Equal Time Correlation Function 



The scaling function for the equal time correlation function is obtained from 
Eq. (H). Hence, 



where 



G{x,t) 



1 



{2ttY Jn Jo 



on on 

1^1 I 1 



dk k'^-^ [i_e-2(s/ifc2+?/2fc4)] [i_cos(A:a)] - 



(32) 



yik'^ + y2k^ 

J33) 

and j = 1, 2. The notation /q means angular integration, and a = cos{k, x). 
We obtain the following power law behavior in the different regimes: 

If pi » p2 then 



G{x, t) oc < 



If pi <^ p2 then 







for X = 0{L) and 






t > ts; 






for ^* ^ a; ^ L and 






t » 


3;2x2(L/2;) 


2X1 


for X L and 












for ^1 » X » ^* and 






tc<^t<^ ts] 






for < ,^1 < X and 






tc<.t-^ ts; 






for X <C ^2 ^ ^* and 












for X » ^2 and t ^ tc- 






for X = ©(i^) and 






t » ts; 




2X1 ^ 


for X ^ L and t ^ ts] 




)2xi; 


for X ^ ^2 and t <^ tg] 






for X ^ ^2 and t <^ ts- 



G{x,t) oc < 



(34) 

To summarize, the correlation function can exhibit both EW and MWV 
scaling behavior for pi ^ P2, while in the opposite case we find "anomalous" 
scaling behavior 



(35) 



of the kind characteristic of the MWV equation [14| as well as certain non- 
linear models |15, 28, 29 1. The scaling (|35|) is anomalous in the sense that 
G at fixed x has no finite limit for L ^ oo and ^,2^00; this implies the 
appearance of arbitrarily large height gradients and is associated with the 
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fact that the MWV roughness exponent X2 > 1 for d < 2, compare to (p^). 
Note that the increase of G with x in (^) is governed neither by xi nor by 
X2, but by an anomalous roughness exponent X = X2 — Xi = ^■ 



3.3 Saturated Correlation Function 

The scahng function for the saturated correlation function is obtained from 
Eq. (H). Hence, 

on 2D 

Cs{x,t) = — = ^X^X^F^ (y^, ys) , (36) 

where 

(37) 

and j = 1,2. For simphcity, we only give the results in the limits where 
either i = or x = 0. In the first case, the power law behavior of Cs{x, 0) is 
given by: 

If pi S> P2 then 

f forx = 0(L); 

Cs{x, 0) oc < x^^^i , for r <C X < L; 

[ x2x2(L/x)2xi , forx«r«;L. (38) 

If pi <^ p2 then 

„ . / forx = 0(L); 

^s[X,U) oc I ^2X2(^/^)2X1 ^ for X « L. 

The behavior of Cs{0,t) is given by: 
If pi » p2 then 

r L2xi, fort>t,; 
Cs(0, t) oc < t^A ^ for tc < i < ts; 

[ t2/32^ fort^tc. (39) 

If pi <C P2 then 

^2,2, fort«ts. 

As in Sec. 3.2, both EW and MWV behavior are found in the different 
regimes. 



4 Applications 

The goal of the present section is to illustrate how linear growth equa- 
tions can be used to address experimentally relevant questions about kinetic 
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roughening that have so far received httle or no theoretical attention. Specif- 
ically, we discuss the contribution of the substrate roughness to the width 
of a growing surface, and the detailed form of the spatial height correlation 
function. 



4.1 Effects of Substrate Roughness 

In the real world, thin films are rarely deposited onto a perfect, atomi- 
cally flat substrate. Consequently every experimental investigation of kinetic 
roughening has to deal with the substrate contribution to the roughness of 
the film surface. To the extent that this problem has been addressed at 
all, it is usually assumed |3C] that the substrate gives an additive, constant 



contribution to the variance of the height fluctuations (the square of the 
width), as 

w\t)=wf + wUt) (40) 

where Wi is the width of the substrate surface and wq denotes the growth 
induced contribution (throughout this section the lateral system size is taken 
to be infinite and therefore the dependence on L is suppressed). This simple 
ansatz ignores the fact that the memory of the initial roughness is lost during 
the growth process, as the growing film successively covers up the features of 
the substrate. The short wavelength features are preferentially suppressed 
(compare to a layer of snow covering a rugged landscape), an effect of much 
importance for the scattering from multilayer films []3l| . 
In the following we show that (i) within the framework of linear growth 
equations, the superposition ansatz (^) is justified, however (ii) the sub- 
strate contribution Wi becomes time dependent and decreases with i in a 
manner governed by the ratio of the substrate correlation length to the 
correlation length ^(t) ~ t^/^ of the growth process; for long times, ^ ^ ^s, 
we find 

u;,(i)~w;,(0)(6/e)"/'- (41) 

The fact that the total width is the sum of a decreasing and an increasing 
part entails the somewhat counterintuitive possibility that w{t) may ini- 
tially decrease with increasing film thickness, as has been observed in recent 



experiments ]3^, |33[ . 

To justify statements (i) and (ii), we assume that the surface has been grown 
from time t = —to to t = driven by some initial noise ?7_ in such a way 
that the height-height correlation function is given by 

{Hoik)Hoik'))^_=fik,k'), (42) 

where Ho{k) = HQ{k,t = 0), and the average has been taken over all con- 
figurations created by the noise rj-. We note that for t > 0, r/_ = 0. At 
t = 0, a new growth process with noise rj is turned on. Let us denote the full 
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solution of the CLG equation by H{k,t) with the initial condition HQ{k). 
Then, 

H{k, t) = i?o(^)e""(^^* + h{k, t) , (43) 

where a{k) = vik"^ + |z^2|^^) and h{k,t) denotes the solution of the CLG 
equation with a flat initial condition. Calculating the correlation function 
with respect to both the old noise t?- (from t = —to to t = 0) and a new 
noise r] (for t > 0) gives 

{{H{k,t)H{k',t))f,)f,_ = 

f{k, k')e'^-('^^+-('^'^^' + {h{k', t)h{k, t))f, , (44) 

given that either (i?o(^))r)_ = 0, or {h(k, t))^j = 0, which makes the crossterms 
disappear. The influence of the initial conditions vanishes exponentially fast 
in the Fourier space but not necessarily in the real space. To see this we 
consider the following example. At time t = —Iq the surface is flat. Then 
we switch on the beam and let the surface evolve until t = driven by the 
CLG growth dynamics (the noise r/_ is white). For / we get 

f{k,k') = {f dTfi^{k,T)e''^^^^ f dT'f}.{k',T')e<^^^')f,_ 
J —to J— to 



2D 1 - e-2«{'=)to 



6'^{k + k') . (45) 



{^T^Y 2a{k) 

Integrating the first term on the right hand side of (p^ ) with respect to k 
shows that the initial roughness decays slowly, as wf ~ to/i, in the limit 
i — > oo, when the substrate dimension d = 2. This corresponds precisely to 
(13) with EW scaling, ~ and ^ ~ t^/'^. 

In the preceding we dealt with the full CLG equation. Next we will focus on 
calculating the time dependent substrate surface width wf (t) for the special 
cases 1/2 = (EW equation, z = 2) and z^i = (MWV equation, z = A) 
of @). Moreover, we assume that the substrate surface has been grown by 
either EW or MWV dynamics, so that its correlations can be described by 
the Fourier amplitudes 

{\Ho{k)W. = ^[l - e'^'^^n (46) 

where is the substrate correlation length, and = 2 or = 4 for surfaces 
generated by EW or MWV dynamics, respectively. The roughness exponent 
of the substrate is Xs = {zs—d)/2. Solving (|3) subject to the initial condition 
Ho{k) one obtains 

{mk,t)\\h^ = {\Homv-e-'''''^^' + ^^[1 - e-^'^^Wn (47) 
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where z is the dynamic exponent of the growth process and ^(t) = (2z^f)^/^ 
its correlation length. Since the variance is obtained by integrating (^) 
over fc, the decomposition (^0|) is valid and the substrate contribution at 
time t is given by 

wl{t) = l^d-i^o dk k'^-^-'^ [1 - e-^'^^^)"] e-^^^^^ . (48) 

Jo 

The upper cutoff A^max is of the order of the inverse lattice constant; in the 
following we assume that and ^ are large compared to l/^max and set 
fcmax = oo in (H). 

Substituting q = (k^sY' in (^) yields 

/■oo , 

w^t) = (l^rf_i^o/^s)e'' / dq q^/^^-\l - e-'^)e-^'^'''' (49) 



where the dimensionless time variable r = has been introduced (re- 

call that ^ ~ t^/^). The integral can be explicitly computed if z = Zg', details 



are given in Appendix A. To analyze (49) in the general case, let us first 
assume that Zs > d, i.e. Xs > 0, which covers MWV substrates in d = 1 and 
d = 2, and EW substrates in d = 1. 

This ensures that the integral over g<^/^s-2 converges at infinity, and hence 
the factor exp(— rg^/^') can be dropped when r ^ 1. The width is therefore 
independent of r for r <C 1. Physically this simply reflects the fact that, 
for a surface with power law roughness, Xs > 0, the width is dominated by 
the long wavelength fluctuations with wavelengths of the order of (,s- At 
time t initial fluctuations of wavelengths up to the correlation length ^(t) 
have been eliminated, thus the substrate contribution to the width decreases 
appreciably only when ^ ~ or r ~ 1. 

For Zg = d (i.e. xo = 0, for example an EW substrate in d = 2) the integral 
over q'^/^n-'i diverges logarithmically at large q. The factor exp(— rg^/^") 
then has to be retained, and one finds that Wi{T)'^ ~ ln(l/r) for r ^ 1; 
of course this behavior is valid only for r > (fcmax^s) ^) since the initial 
roughness tfj(O)^ f« Q-d-iA^ ln(A;niaxCs)- Finally, in the (somewhat academic) 
case Zg < d, Wi{T) decreases as a power law also for r ^ 1, as Wj(r)^ ~ 
^-(d-zs)^ see Appendix A. 

The behavior for large r can be discussed independent of Zg- For r » 1 the 
integral (49) is effectively cut off at g ~ -y-"^"/^ <C 1, therefore we can set 



1 — e ~ g in the integrand, and it follows by rescaling that 

wfir) ^ Aoe~'(6/0' ~ t-''/\ (50) 

This result has a simple interpretation. For ^ ^ the growth process av- 
erages over a large number N = (^/^s)'^ of domains in which the initial fluc- 
tuations are statistically independent. The height variance of each domain 
is of the order AQ^g^~'^, and averaging over domains reduces the variance 
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by a factor of If the substrate has power law roughness {xs > 0), then 
Wi{Of « Ao^l--'^ and (|5^) is identical to (0)- 

The last argument is valid also for completely general initial conditions char- 
acterized by a correlation function of the form 

(\Ho{k)\\. = "^ms) (51) 

where xo denotes the roughness exponent of the substrate and the scaling 
function g{s) satisfies g{s —>■ oo) = 1, g{s — > 0) = 0. The expression (^) 
generalizes to 

/■oo 

wfir) = 17,„iAoe,'» / dq q'^'+^^^^g{q)e--'^\ (52) 

To extract the behavior for large r we need to know how g{s) vanishes for 
small s. This is fixed by requiring that (51) should have a finite limit for 
A: — >■ (this limit gives rise to the center of mass fluctuations of the surface, 
see iQ). Consequently g{s) ~ s^^+^xo^ ^nd (|^) decays as 

wKr) « n^-iA^il^H-'"^ « wKms/if (53) 
in accordance with the heuristic argument. 

In summary, we have shown that, under rather general conditions, the sub- 
strate contribution Wi to the width of a rough growing surface remains es- 
sentially constant as long as the correlation length ^ of the growth process 
is smaller than the substrate correlation length ^g-, and that it decreases 
according to ( |4l| ) for ^ ^ ^s. While the functional form of the transition 
between the two regimes is not analytically accessible in general, a useful 
interpolation formula is 

Wi{tf = Wi{<df{l + t/ts)-'"' (54) 

with a fit parameter t^, which has already been employed in the analysis 
of experimental data |^^. In Fig. 2 the decay of the substrate width is 
illustrated for two special cases, and the formula ( |5^ ) is compared to the 
exact expression derived in Appendix A. 



4.2 Shape of the Height Correlation Function 

The dynamic scaling hypothesis of kinetic roughening theory [^l| states that 
the height -height correlation function should have the scaling form 

C(x, t) = {h{x + X, t)h{x, t)) = w'^{t)F{x/C{t)) (55) 

with w"^ = C{0,t), hence F(0) = 1, and ^ ~ t^/\ To date almost aU theo- 
retical work has focused on the behavior of the height difference correlation 
function 

G{x,t) = {ih{x' + x,t) - h{x' ,t)f) 

= 2[C{0,t)-Cix,t)]=2w^[l-F{x/m)] (56) 
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for X <C ^, i.e. the short distance behavior 

F(s) ^ 1 - 0(s2x) ^ s^O , (57) 

of the scahng function. In contrast, the overah shape of F or, in particular, 
the way it decays for large arguments has not been addressed theoreticahy, 
although a considerable amount of empirical information is available pO|| . 
This question is of some experimental importance, since the interpretation 
of scattering data from rough surfaces typically requires a model for the 
entire correlation function |3^. A widely used form for F is [^] 

F{s) = exp(-s2>^) , (58) 

which assumes a simple relation between large distance and short distance 
behaviors. Here we show by explicit calculations within the linear theory 
that the large distance decay of F involves a new exponent not simply related 
to and that moreover the decay can be oscillatory rather than monotonic. 
First we make some general remarks concerning the correlation function for 
the CLG model. The entire correlation function valid for all ranges can be 
rewritten for d = 2 as (by substitution k' = k/^ in Eq. (^)) 

,2 ,4 

G'(f'|.|) = — I / dkk[l-Mkx/OV--, J, (59) 

where = {li^it^^'^ and ^2 = (2|i^2|t)^^'^. As far as the CLG model can 
be considered to describe the physics correctly either Eq. ( |59| ) or Eq. (^) 
could be used directly as fitting functions, with the parameters {D, z/i, 1/2) or 
{A = £)/(27rz/i), ^1, ^2) to be fitted. The amplitude is set by parameters D or 
A. The correlation length ^(t) in Eq. ( [59| ) corresponds to the experimentally 
observable correlation length. We can numerically determine the dependence 
of ^ on ^1 and ^2 by using the relation C{x = ^)/C{x = 0) = e'^, which is 
the usual definition of the correlation length. The surface ^(^1, C2) is plotted 
in Fig. 3. A least squarest fit for a set of solutions ^ = (,{Ci,(,2) yields 

= ai6 + 026 = 0.085336 + 1.605736 • (60) 

In addition, numerical studies show that ^ follows the scaling relation ^(a^, 06) 
a^(6j6)- Our choice (,p = ai6 + 0^26 is a simple case which quite accu- 
rately satisfies the definition for the correlation length. 
Unlike in the case of the heuristic correlation function (|58[), the requirement 
X <C ^ alone is not sufficient to determine the scaling behavior of G for CLG 
model with a fixed exponent x- The scaling analysis of Sec. ^ shows that the 
scaling of the correlation function depends on two dimensionless quantites: 
6/6 ~ (t/tc)^^^ and ^ix/^^ ~ x/e*. If t/te > 1, x/r > 1, and x < 6 
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then G oc x2xi('^=2) and x ~ Xi{d = 2) = 0. If t/tc <. 1, x/t <C 1, and 
X <C then G oc x2>^2(c(=2) ln(^2/x) and thus x ~ X2(d = 2) = 1. 
Let us next determine the behavior of the correlation function for large 
values of the ratio x/^ = s. Rather than trying to present the full solution 
for the CLG model we shall limit ourselves to the one and two dimensional 
EW and MWV models. For the d = 1 EW case, the evaluation of C is 
straightforward. We have 

dtC[x,t) = ^ dke'^'-^ . (61) 



This is a Gaussian integral which can be done. The integration over t is also 
standard. We obtain the result: 

C{x,t) =w'^{t)F{s) , (62) 

with 

w\t) = (63) 

F{s) = - [' dyy-y\~^'li^y) ^2s~^e-'"l^ (64) 
2 Jo 

as s ^ oo. The behavior of the scaling function F is in striking contra- 
diction to the conventional ansatz (^), which would lead one to expect an 
exponential decay for x = 1/2. 
The d = 2 EW case is also simple to calculate: 

dtC{x,t) = — r dkkJo{kx)e-^^^^'^^ . (65) 



TT Jo 

Integrating the result from to t yields |13] 



n /-oo T-) 



^(^) = -A / ^ 1 ^1(^74) , (66) 

where Ei is an exponential integral. The asymptotic power law expansion 
of El reveals that as a function of s the leading term is equivalent to the 

d=l case: Ei{x) ~ x~^e~^'(l — x~^ + 2x~^ ), hence F{s) ~ s~^e~* 

In the d = 1 MWV case we confront the integral [p^ j 

/oo , /-oo 

dke''^-^ =i:^^ dke-f^^K (67) 
-00 J ~oo 

The behavior of the leading asymptotic behavior can be evaluated using the 
method of stationary phase [^] . In the limit s — > 00 only the neighborhoods 
of saddle points contribute to the integral. They are determined from the 
relation dkf = 0: 
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Removal of each saddle point to the origin results in a quadratic dependence 
of f{k) on k: f{k) — f{ki) bik"^ as ^ 0, i = 0, 1, 2. Thus, in the leading 
approximation we retain only the Gaussian part of the integrand: 

dtC{x,t) oc i2\^'^(i'^^^'^ ■ (69) 

Furthermore, we note that the contributions coming from distinct saddle 
points are additive if the paths of integration are properly chosen in the 
complex plane. Despite that the individual contributions coming from points 
ko and ki are complex, their sum is real. The point ^2 has to be discarded 
as it would cause an exponential blow-up. Finally, we obtain 

C{x, t) ~ ^s"5/3 g-ai«*/3 gij^(^2^4/3) ^ ^2 cos(a2S^/3)] , (70) 

where oi = (3/8)4^-'^/'^ and 02 = (3\/3/8)4^^/^. The correlation function C 
is an oscillating function for large values of the argument s and can assume 
negative values as well. 

For the d = 2 MWV case we use the same method as in the previous case: 

o n roo POO 

d^Cix,t) = -0^ J dkij dk2e-^^2ikI+4fMkix^+k2X2) ^ (71) 

To determine the saddle points we set d^^f = dk^f = 0, where / denotes 
the expression in the exponential of ([7l|). This yields 

dkj = ; dkj = 0, (72) 
AClkf + A^lkikj -ixi = ; 4C|A;| + 4^|A;2A:? - 1x2 = , (73) 

which can be solved with the ansatz ki = axi, k2 = ax2- Again we get 
three saddle points, one of which corresponds to an exponential divergence. 
Summing up the contributions of the remaining two and integrating over 
time gives a real answer: 

C(s) ~ ii- e-'^^^'" cos{a2s'/') . (74) 

As before, an oscillating function results. 

Summarizing, it appears that the large distance decay of the scaling function 
F{s) in ( |55[ ) generally has the form of a 'squeezed' exponential decorated by 
a subleading power law, 

F{s) ^ s-^'e-"'', s^oo, (75) 

where the amplitude c can be real (as in the EW equation) or complex 
(as in the MWV equation); in the latter case the decay is modulated by 
oscillations. In contrast to what is suggested by the heuristic function (^), 
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the decay exponent has no direct relation to the roughness exponent x\ 
rather, in the examples treated here, 6 seems to be characteristic of the 
surface relaxation process but independent of dimensionality. In Appendix B 
we give an alternative derivation for d = \ which suggests that 5 = z/{z — l) 
and 7 = 1 + 5/2 for the generalized linear equation 

dth = -{-dly/'^h + r] (76) 

with integer values of z/2. 

It is worth pointing out that a rapid, exponential-like decay of the correla- 
tion function should be expected only for growth processes that are local in 
space, such as (^), for which the Fourier transform of C{x,t) is an analytic 
function of k for small \k\. A simple example of a nonlocal kinetic rough- 
ening process is diffusion limited erosion |3^], formally the generalization of 
(|7^ ) to z = 1, for which it can be shown that in d = 1 

C{x,t) = ^-lnil + e/x') (77) 

for distances larger than the lattice spacing a; at x ~ a, C ~ ~ 
(L>/i^)ln(^/a), since x = (1 - d)/2 = in d = 1 [||. As usual, the cor- 
relation length is ^ = 2i't. At large distances ( [77| ) decays as a power law, 
C ~ (i^/x)^, reflecting the nonlocal nature of the interface dynamics. 



4.3 Effects of Noise 

We have shown in Sec. § that the correlation function for the noise typically 
consists of two components, 

(r?(f , t)ri{x, t')) = 2{Di - D2V^)5{x - x')S{t - t') (78) 

where the first, white noise component arises from evaporation and deposi- 
tion, while the second, conserved component reflects the thermal fluctuations 
in the surface current pi. In Fourier space the correlator (|78D reads 



(r?(^, uj)fi{k\ J, ))f, = {2t:)-\2Di + 2D2k'^)5\k + k')5{uj + J) , (79) 

where r) denotes the Fourier transform with respect to both space and time. 
So far we have considered the behavior of the surface on scales large com- 
pared to the noise crossover scale I** ~ {D2/ Di)^/"^ , where the conserved 
part can be neglected. In the opposite regime x <C I** we can set Z)i = 
and analyze the behavior generated only by diffusion noise. The correspond- 
ing power laws are easily evaluated: Because the D2 correlator in ([79|) is of 
the same form as for the white noise apart from the factor k"^, the results 
presented in Sec. |3| hold if we replace dhy d + 2. 



In Ref. 1 41] it was pointed out that the effect of fast degrees of freedom may 



sometimes be represented by means of a coloured noise term r/c, i.e. the noise 
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has power law correlations. Through renormalization group analysis it was 
shown that power law correlated noise changes the scaling exponents for the 
KPZ equation. Now we shall extend our analysis of Eq. (jl]) to include power 
law correlated noise as well. We assume that the noise has long distance 
correlations of the form 

{rjc{x,t)vc{x',t'))^^ = Ix-xfP-^lt - tf^-' , (80) 

with < p < d/2 and < < 1/2. Taking the Fourier transform of Eq. (^) 
leads to 

iUk, a;)r/(fc', u;^^ = ^2^6''{k + + u') , (81) 

where D{k,tLi) = D'\k\~'^P\Lo\~'^^ , and D' is a function of p and g. The 
analysis of the power law behavior for both EW and MWV cases we get the 
following relations between the new roughening exponents (primed) and the 
old ones: 

f3[ = Pi + lp + Q ; x'i = xi+P + '^e ■ (82) 
P'2 = P2 + ^P + Q ; X2 = X2 + p + ^g- (83) 

Clearly, z[ ^ x'l/P'i =2 and z'2 = x'2/ — satisfied because Xi/Pi =2 



and X2/P2 = 4. These results agree with previous treatments |12, 16, |4l[|. 



5 Summary and Conclusions 

In this work, we have presented a detailed quantitative analysis of a com- 
bined linear stochastic growth equation, which incorporates both EW and 
MWV type of behavior. In particular, we have identified the relevant scaling 
variables and calculated in detail the behavior of the surface width, and var- 
ious correlation functions in different regimes. Our analysis shows that the 
CLG model possesses MWV type of scaling regime at early times {t <C tc) 
and crosses over to EW type of growth for later times. In a finite system 
it is, however, possible that the system maintains MWV type of growth 
dynamics for all times given that the prefactor of the fourth order gradient 
term is very large (L <C £*)■ Additionally, the equal time and saturated 
correlation funtions reveal interesting and complicated crossover behavior 
between the two limits. We should also note that in the appropriate limits, 
our results completely agree with the previously derived results for concern- 
ing the surface width and the equal time correlation functions of the EW 
and MWV equations |, |l|, |l|, |T|, |l|, |l^. We have also recalculated the 



scaling exponents in Sec. ^ for coloured noise that has long range spatial 
and temporal correlations. Our results show that the long-range part of the 
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noise changes the exponents, in analogy to the case of the nonhnear KPZ 
equation plf . 

The solvabiUty of the CLG model allows us to address other important prob- 
lems related to surface growth problems as well. A particularly interesting 
case concerns the influence of rough initial conditions. In Sec. ^ we demon- 
strate that the effect of initial roughness in d = 2 disappears as a power 
law in time ~ to/t if the surface has initially roughened through the CLG 
growth dynamics for a time to- In general, the initial roughness vanishes as 
i^s/O'^ ~ for EW and MWV models as ^ ^g. For a general growth 

process the behavior of this transient is easy to evaluate by Fourier trans- 
forming the exponentially damped initial height-height correlation function 
(Eq. (0)). 

Finally, to make the connection to experiments more concrete we have pro- 
posed the use of the parametrized correlation functions obtained for the 
CLG model for fitting of diffusive x-ray refletivity data. In Sec. ^ we estab- 
lish a relation between the parameters vi and U2 of the CLG growth model 
and the experimental fitting parameter ^ (cf. Eq. (|60[)). Moreover, both 
the small [x <C and the large (x ^ ^) scale behavior of the correlation 
function of the CLG model are found to be different from the experimentally 
used fitting functions. This also implies that the measured decay exponent 
X from heuristic fitting functions bear no direct relation to the roughness 
exponent of the CLG model. We hope that these calculations will be useful 
in future experiments on surface growth. 
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Appendix A: Substrate Roughness for Zg = z 



Here we explicitly evaluate the integral (|49| ) in the special case Zg — z. 
Taking the derivative of ( ^9[) with respect to r and reintegrating we obtain 



„2 _ ^d-1^0r(d/ z) cz-du I 1 \l-d/z _ ^1-d/z] 



d 



-is[{T + ly-"/' - T^-^i^i (84) 



Three cases have to be distinguished. 

(i) z > d: In this case both the substrate and the growing surface are rough 
with roughness exponent x = (z — (i)/2 > 0. The expression ( |8^ then has 
a finite limit w^fS) = 17rf_iAor(d/z)^|~'^/(z - d) for r ^ 0. For r <Cl the 
roughness remains essentially unchanged, 

w.irf^wmH'^-r'-'"'), (85) 
while for r ^ 1 it decays as 

w^,{T)^{l-d/z)wKQ)T-'/^ (86) 
in accordance with (^T]). 

(a) z = d: Substrate and surface are logarithmically rough. Taking the limit 
(z - d) ^ in (H) gives 

wUr) = (nd-iAo/z) ln(l + 1/r) (87) 

which diverges for r — > 0. Clearly (^) is valid only for ^ > i.e. 
r > To = (/cmax^)"''; the initial roughness is Wi^To)'^ « Qd_i^o ln(fcmax6)- 
For tq <^ T <^ 1 the substrate contribution decays slowly, as Wi{T)'^ 
f^d-i^o lii(l/r), while for r » 1 ( |87| ) reduces to 

-1 



u;f(r) « (l^d_iAo/z)r 

which is of the form (^) with a prefactor of order unity. 

(in) z < d: Substrate and surface are smooth, and the width is microscopic 

at all times. In this case (jsj) combines two power law decays: For r ^ 1 



independent of ^s, while for r ^ 1 one finds 

wKr) « [17d_i.4or(d/z)/z]C('-^V-^/^ ~ if^-'Kis/O" (90) 
as in the previous two cases. 
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Appendix B: Derivation of the Exponent S 



Here we present a simple way of computing the exponents 5 and 7 char- 
acterizing the decay of the height correlation scaling function (^5|) for the 
generalized linear equation (|7^) in one dimension. It follows immediately 
from d?^) that the Fourier transform F{q) of F{s) is, up to some factor, 



Fiq)^q-^{l-e-'''). (91) 
We now exploit the fact that by definition, F{q) can also be written as (Eol 



00 . 

m^T.^i-'iT (92) 



n=0 



where 

/oo 
ds s^'Fis) (93) 
-00 

is the n^^ moment of the real space scaling function. Comparing ( |9l| ) to the 
power series (p^), we obtain the relation 

, . , 1 {zkV. , , , 

The asymptotics of F{s) can then be extracted by comparing the behavior 
of (Q) for large k with that obtained from the ansatz ( [75|) inserted into 
(j9^); assuming the amplitude c in ( [75| ) to be real, the integral ( p3[ ) is easily 
evaluated in the saddle point approximation. The leading behavior is An ~ 
77,"/'^ which, compared to \Azk\ ~ }^{^~'^)^ from ([9^), implies that 5 = 2;/(z — 
1). Moreover, from the integral (^) one obtains a power law factor ~ 
j7,i/'5-i/2-7/'5 which has to be matched to the factor (k + l)""*^ that appears 
in (|9^ ; this forces 7 to take the value 7 = 1 + 5/2. Both results agree with 



the detailed calculations of Sec. 4.2 
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Figure Captions 



Fig. 1 A log-log plot of the surface width w'^{t) of the CLG equation for 
several values |z/2| and z^i in d = 1, with L = 10000. The two scaling 
regimes of the surface width w'^ oc 

for t <C tc and 

for t ^ tc are clearly visible. For the topmost curve ui = 1 and 
|z^2| = 10~^ and the crossover time Intc = ln(|i/2|/i^i) ~ — 18 is shown. 
For the second curve z/i = |i/2| = 1, and Intc = 0. Finally, the two 
overlapping lowest curves represent the CLG surface width with vi = \ 
and 11^21 = 10^, and the pure MWV case with vi =0, |i^2| = 10^. In 
both cases the saturation time Ini^ ~ 18. 

Fig. 2 Decay of the substrate contribution to the surface width for dimen- 
sionality d = 2 and dynamic exponents Zg = z = A (full upper curve) 
and Zs = z = 2 (full lower curve). The full curves show the exact 
expressions derived in Appendix B, while the dotted curve shows the 



interpolation formula (54). The dashed lines indicate the asymptotic 
power laws. 

Fig. 3 (a) The correlation length ^ as a function of the two dynamical correla- 
tion lengths ^1 = (2i/it)V2 and 6 = {'^\^2\tY'^ of the CLG model, (b) 
The relative error E = (^p — ^)/^ between the plane fit = 01.^1-1-02^2 
and ^ (Eq. ([5^). The average relative error is about 1%. The maxi- 
mum error in the figure (10%) occurs for large ^1 and small ^2- 
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Figure 2: (Majaniemi, Ala-Nissila, Krug) 
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Figure 3: (a) (Majaniemi, Ala-Nissila, Krug) 
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Figure 3: (b) (Majaniemi, Ala-Nissila, Krug) 
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